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This paper derives a slope overload noise power formula for linear 
delta modulators having ideal integrators and Gaussian random inputs. 
Although the same problem has been treated by others, the present result 
is the only one applicable to all slope-following capacities and input spectra. 

Despite its singleness of purpose, the paper divides logically into two 
parts. In Part 1, a common element in all previously published results 
is used to derive a new slope overload noise power formula. This derivation 
is analytically rigorous and provides some useful insights, but pertains 
to a particular kind of spectrum and so is incomplete. 

The more universal result we seek is derived in Part 2. The approach 
here is far less rigorous and amounts to approximating the influences of 
other kinds of spectra by modifying the result of Part 1 . The final expres- 
sion contains four spectrum-related coefficients, for which simple formulas 
are given, and has an estimated accuracy of 1 dB for all cases of practical 
interest. Computed results are given for two important families of spectra 
and comparisons are made with previously published results. 

Part 1 

I. INTRODUCTION 

l.l Objective 

This paper presents some new theoretical results on slope overload 
noise in linear delta modulators. In particular, we derive the mean 
power of the (unfiltered) slope overload noise at the demodulator 
output when the modulator input is a stationary Gaussian random 
process and the modulator feedback path contains an ideal integrator, 
Fig. 1. This same problem has already been treated by other investi- 
gators, notably Zetterberg, 1 Rice (with O'Neal), 2 Protonotarios, 3 and 
Abate,* but several factors have prompted a reexamination. One is 
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Fig. 1— Linear AM codec with perfect feedback integrator, (a) Equivalent block 
diagram, (b) Input and feedback signals, (c) Model used. 

that some important disparities among the results of these separate 
studies remain unresolved ; another is that the existing formulas do not, 
either individually or collectively, pertain to all slope-following 
capacities and input spectra ; and finally, few explicit clues are avail- 
able as to just how accurate each formula is and where it loses validity. 
In response to this situation, we have endeavored to find an expres- 
sion for slope overload noise power that is accurate for all slope- 
following capacities and input spectra of possible interest. The result 
reported here satisfies that objective. 

1.2 Noise Descriptions and Definitions 

The idealized AM codec (coder/decoder) shown in Fig. 1 ex- 
emplifies the process we want to analyze: Every t seconds the in- 
put signal (x(t)) is compared with a locally quantized version of 
itself (y(t)), and a unit impulse is generated with a polarity that is 

/positive \ . f (*>V\ The resu iting binary impulse stream is 
^negative/ \x < y) 

applied through a feedback gain factor (5) and ideal integrator to 
produce y(t). At the decoder, y(t) is reconstructed by using a replica 
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of the coder feedback network, and a final low-pass filter smooths the 
sharp edges of y(t), yielding a closer approximation to x(t). 

Since y(t) cannot change by more than 5 units in t seconds, 5/t is 
the highest input signal rate-of-change that the AM codec can follow. 
We call 5/t the AM slope-folloicing capacity and denote it by x' . When 
\dx/dt\ exceeds this quantity, slope overload occurs and gives rise to 
the kind of error shown in Fig. lb. In addition to this sporadic form 
of distortion, the "hunting" of y(t) for x(t) by means of quantum steps 
gives rise to a perpetual distortion called granular noise. Obviously, 
granular noise is reduced by decreasing 5, but at the expense of a 
reduced slope-following capacity and, hence, greater slope overload 
noise. 

To delineate slope overload and granular noise for purposes of this 
analysis, let us suppose that x(t) is passed first through a AM codec 
having infinitesimal 5 and t, but with the ratio between them the same 
as in the actual system (Fig. la). The decoded output (before the 
final filter) will then be the smooth function y(t) shown in Fig. lc. 
If y(t) is then passed through the actual system, the decoded signal 
will be a very close approximation to y(t), Fig. lb. In agreement with 
previous conventions, we define slope overload noise to be the differ- 
ence between x(i) and y(t) ; and the remaining distortion, (y(t) — y(t)), 
is defined to be granular noise. This essentially is the approach tacitly 
followed in most of the published literature on AM noise. 1-7 

We define a slope overload noise burst to be the nonzero difference 
between x(t) and y(t) over an interval [t a ,t b ~], which starts because 
\dx/dt\ > x' &t t = t a and ends because y(t) intersects x(t) at t — t b , 
(e.g., see Fig. 2, in which t a = 0). The mean power of these bursts, 
averaged over all time, is the slope overload noise power. 

Two qualitatively different kinds of noise bursts can be identified. 
One is the kind initiated when dx/dt\ increases through x' while y(t) 
is following x(t), which we call primary noise (see Figs, lc and 4b). 
The other arises when a prior burst terminates at a point where | dx/dt \ 
already exceeds x'„ (see Fig. 4c) ; the new burst that commences at this 
point we call secondary noise. 

Finally, the slope overload factor is defined to be the ratio of the 
slope-following capacity to the rms input signal derivative, i.e., 

A Xo 

S = (1) 

(dx/dt) v ' 

rms 

It should be obvious that the noise power decreases monotonically 
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Fig. 2 — Zetterberg's noise model. 



with increasing S. It is also clear that when S is large (in which case 
\dx/dt\ rarely exceeds x' ), secondary bursts are rare so that primary 
noise dominates the noise process; and that, by the same token, 
secondary noise dominates the noise process when S is small. 

1.3 Equivalent Input Process 

The generality of our analysis will be enhanced if we can assume 
that |.r(0l is a bandlimited process. This assumption is made valid 
by regarding ( x (t) } as an equivalent process related to the true one as 
follows : Let the true input process be { x (t) } , having a power spectrum 
X (f). Since the delta modulator really acts on discrete samples of 
the input separated by r seconds, its response is the same as if the 
input were a process { x (t) } having sample functions of the form 



00 sin (ir(t — fcr)/r) 

x(t) = £ x (kr) 

*=-» (v(t - kr)/r) 



(2) 



and a power spectrum 



Vi/i - E A (/ + -); ^ /£— only. (3) 
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We will assume here that the sample function and slope overload noise 
depicted in Fig. lc correspond to the bandlimited process {x(t)\, i.e., 
that the input spectrum is X(f). If \x (t)} is bandlimited by some 
frequency W ^ 1/(2t), the two processes, and consequently their 
spectra, are one and the same. (This condition is usually tacitly- 
though not always rightly-assumed in AM noise analyses.) If (rr o (0| 
is not so limited, the difference between the two processes is aliasing 
distortion, which can be analyzed separately. 

1.4 The Spectral Moments 

The spectral moments of the input process play a decisive role in 
establishing the past and present slope overload noise formulas. 
Following the convention of most authors, we define the nth moment 
to be 



b n = I u 2 »X(f)df. (4) 

The bandlimited nature of X(f), as discussed above, guarantees the 
finiteness of all the moments. It is easy to show that the complete 
set of these finite moments, (b , b h b 3 , ■ ■ •, b„, • ■ •), uniquely deter- 
mines X(f), just as the converse is true. 

It should be obvious that b„ is the mean power of the nth derivative 
of x(t), i.e., b n = (x {n) (t)) 2 ': Thus, for example, V&7 is the rms derivative 
of the input and so (1) can be rewritten as 

S = x' /^b~ 1 . (5) 

Also, bo is nothing other than the input signal power, as distinct from 
the ac input power which we denote by a 2 . If the input contains no 
discrete dc component, then 6 = a 2 ; otherwise, b exceeds a 2 by the 
amount of the dc power. 

Finally, we note from (4) that b n is real and non-negative because 
X(f) is. An additional relationship for b„, derived by applying the 
Schwarz inequality to (4), is 

b n ^ — ; n ^ 2, (6) 

0„_2 

that is, b n for n ^ 2 has a lower bound determined by the previous 
two moment values. We shall utilize this relationship later. 
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1.5 Outline of the Work 

Our objective is to find a suitable approximation to the true slope 
overload noise power formula, the latter being denoted by N(S). 
Section II describes the methods and results of previously published 
analyses. While the various approaches differ, they all suggest that 
this noise power is essentially determined by S and the first three 
spectral moments, b , b lf and b 2 . Section III treats this possibility as a 
premise and identifies a simple two-band process, the parameters of 
which can be chosen to yield any (6 , b u b 2 ) and which can be analyzed 
precisely. The slope overload noise power for this process is derived 
and is denoted by N t b(S). 

Unfortunately, N tb (S) is not precisely applicable to all spectra 
having the same b , b h and b 2 . Instead, it is a lower-bound variation 
for all such spectra, which becomes increasingly unreliable in general 
as S decreases. To obtain a more universal estimate [denoted by 
N (S)2 it is necessary to determine how much the slope overload noise 
power deviates from N, b (S) for spectra other than the two-band kind. 

To this end, Section IV derives a noise power formula, N t (S), 
applicable to all spectra but confined to the large-£ region (S ^ 3.5), 
and Section V derives a noise power formula, N,(S), applicable to all 
spectra but confined to the small-/S region (S « 0). Section VI then 
combines these results to obtain a two-region approximation, N„(S), 
which is accurate for all X(f) and S. The result is given by (62) and 
(77), where a h a 2 , a h a 2) a 3 , and a 4 are related to the spectrum parame- 
ters bo, h, b 2 , b 3 , b<, 6 5 , X(0), and />- 2 [X(/) - X(0)]d/. Further 
study shows that (77) alone can be used over all S with an accuracy 
of 1 dB or better up to S = 6.5, beyond which point N(S) is at least 
119 dB below the input signal power. For practical purposes, therefore, 
our final approximation to N(S) is just (77), with Oi, a 2 , a 3 , and_a 4 
given by (72) through (75). 

Section VII demonstrates the new result for two important families 
of spectra, and compares N„(S) with the noise power formula of 
Protonotarios. The latter is found to be highly accurate for most 
spectra of practical interest and S ^ 2.0. 

II. REVIEW OF PREVIOUS WORK 

In the approaches of Zetterberg, Rice, and Protonotarios, #(&) is 
approximated as the product of the mean energy (8) per slope over- 
load burst and the average rate of occurrence (R) of such bursts. A 
burst is assumed to commence whenever |x'(0| increases through 
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the value x' 0> which means that only primary noise is considered in 
these studies. The average rate of these events is known from the 
earlier work of Rice 8 to be 

R = -J b ^exp{ -(x'.)V26i}. (7) 

IT * Ol 

To find 8, Zetterberg uses the model shown in Fig. 2, where the solid 
line of slope x' represents the demodulator output for the duration of 
the burst, i.e., from t = to t = tb. Clearly, the noise is 

n(t) = x(t) - 0(0) + x'oQ) ^ t g t b only (8) 

which is the first burst depicted in Fig. 2. To avoid deriving the 
random quantity tb, Zetterberg regards the slope overload burst to be 
the entire excursion of x(t) above the line [.r(0) + x' Q, which is the 
variation m(t) in Fig. 2 rather than just n(t). This is obviously an 
approximation, since m(t) can contain one or more spurious bursts 
that are not really part of the original noise burst, as seen. Such 
bursts, however, are low in both energy and probability of occurrence 
when S» 1. 

Zetterberg proceeds by finding, at each t, the average of m 2 over 
the noise burst ensemble, and then integrating over all time and 
multiplying by R, (7). If done correctly, this leads to an estimate of 
N(S) which is marred only by the inclusion of spurious burst con- 
tributions and the ignoring of secondary noise effects. In fact, however, 
Zetterberg errs in defining the initial conditions of the burst and in 
averaging over these conditions, leading to an incorrect solution.* 
In addition, Zetterberg makes a number of functional approximations 
and at least one algebraic error.* For all these reasons, his final result 
is incorrect and will not be repeated here. 

The approach of Rice reverses that of Zetterberg, in that the energy 
per burst is found first and then averaged over the noise burst ensemble. 
The problem of spurious bursts is thus avoided, but at the expense of 
having to find tb- Rice accomplishes this by expanding x(t) into a power 
series about t = 0, and then ignoring fourth-order powers in t and 
higher, an approximation that gains in validity with increasing S. 
In addition, he approximates the third derivative of x(t) at t = by 
its conditional mean given that x'(0) = x' . This can be shown to be 
— (bz/b )x' o , and the relative variation of .r'"(0) about this value tends 
to be small when S » 1. 



* See Protonotarios 3 for a discussion of these errors. 
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Using these two approximations, and correctly identifying and 
averaging over the initial conditions of the burst, Rice obtains the 
result 

243 ^Uexp(- SV 2). (9) 



N °™ = ^WAVjj> 



Note that the ratio of noise power to signal power (6 ) depends solely 
on S and 7 = bl/b b 2 . From (6) we see that 7^1 and observe that 
7=1 only for an infinitely narrowband spectrum. 

It should also be noted that, as S -> 0, N R (S) increases without 
bound at a rate S~ 5 . This is clearly not a true representation since 
N(S) should approach the ac signal power (o- 2 ) in the limit of zero 
slope-following capacity.* 

Noting the disparity between Zetterberg's result and Rice's, even 
at large S where both should converge to exactness, Protonotarios has 
attempted to settle the issue by combining Zetterberg's more accurate 
model with Rice's correct averaging procedure. His analysis leads to a 
double integral solution which is exact except for the inclusion of 
spurious contributions and the ignoring of secondary noise effects. 
In reducing this formal result, Protonotarios makes some functional 
approximations and obtains the following : 



where 



N P (S) = -^=(£) -^exp (S*/2)A( X ) (10) 

2 5/8 2 1/8 

X = 5/(6i/6 2 bo) (11) 



and A ( • ) is a function involving powers and exponentials of the 
argument, [eq. (66) of Ref. 3]. Once more, the ratio of noise power 
to signal power depends solely on S and 7. In the limit as S — > » , A (x) 
approaches 1 and so N P (S) converges to Rice's result, N R (S). In the 
limit as S — ► 0, A(x) varies as S* so that N P (S) varies as S~ l rather 
than S~ b . This is still an unbounded increase, however, so that (10) 
is still not acceptable as a complete characteristic. 

Abate's derivation of N(S) follows an approach quite different from 
the others. Using simulation results reported by O'Neal 2 for three 
particular spectra, he has developed an empirical relationship between 
the AM sampling frequency and the slope overload factor for which 

* If the input contains a discrete dc component, the delta modulator will follow 
it exactly so long as S ?* 0; hence the result lims-»o N(S) ■ a 1 . 
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granular-plus-slope overload noise power is minimized. Combining 
that relationship with a simplifying approximation to van de Weg's 
formula for granular noise power, 5 Abate obtains the following expres- 
sion for slope overload noise power : 

NAS) = Ifi (2^T (1 + 3S) ^ ( ~ SS) (12) 

where W is the signal truncation bandwidth.* We see that N A (S) goes 
to a finite value as S — ► 0, and that its variation at large S is expo- 
nential; in both respects, it differs from Nr(S) and Np(S). If we now 
define 

a 8t 2 6 2 

K = , (13) 

27 (2irW) 2 b 1 

(12) reduces to the form 

N A (S) = K (jj (1 + 35) exp (-35). (14) 

Evaluating K for the spectra studied by O'Neal, we find that it lies 
between 1.074 and 1.75 for the three cases, a range of just 2 dB. 
It is tempting to speculate, therefore, that (14) is a more correct form 
in general than (12), with K a universal constant on the order of unity, 
and that the apparent 2-dB spread in K over the three spectra is a 
result of experimental uncertainties. If we accept this notion, we again 
have the result that the ratio of noise power to signal power depends 
solely on S and 7. 

III. NOISE POWER IN TERMS OF b , b h bo 

3.1 The Two-Band Process 

The published results cited above suggest that N(S) is essentially 
the same for all processes having the same values for the zeroth, first, 
and second moments. Assuming for the present that this supposition 
is correct, we call attention to the two-band process, the spectrum of 
which is shown in Fig. 3a. In the limit as the two bands become in- 
finitely narrow, the zeroth, first, and second moments become pre- 
cisely 6 , &i, and b 2 . By analyzing this simple process, therefore, we 
can derive an exact noise power formula [A^tOS)] in terms of 60, &i, 
and 6 2 and then see how universal it really is. 

* The three spectra treated by O'Neal are of the truncated Butterworth type, 
with corner frequency-to-bandwidth ratios of 0.068, 0.2o, and 00. The data used to 
derive (12) cover an S-range from about 1.4 to 4.2. 
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Fig. 3— The two-band random process, (a) Power spectrum, (b) Sample function. 



To appreciate the simplicity of the two-band process for purposes 
of analysis, we should view it in the time domain. A sample function 
is shown in Fig. 3b and is seen to consi st of a dc level, D, plus a sinusoid 
having radian frequency u = Vb 2 /bi, phase </> = u t+, and amplitude 
A. D is a Gaussian variate, whose mean and variance across the sample 
function ensemble are and (b — &i/&0, respectively; is a uni- 
formly distributed variate on [-*-, +»]; A is a Rayleigh variate of 
mean-square value 26f/6 2 ; and D, <f>, and 4 are mutually independent. 

We can derive the slope overload noise power for this process by 
finding the mean noise power associated with a given sample function 
and averaging over the distributions on D, <t>, and A. This approach is 
simplified by the fact that D and <f> do not really influence the result. 
That is, the noise pattern for a given sample function converges 
ultimately to a variation about D that depends solely on the sinewave 
amplitude and frequency (Fig. 4). 

We thus see that the noise power per sample function is the steady- 
state noise energy per half-cycle, denoted by S(A), times the rate of 
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Fig. 4 — The three regions of sinewa ve amplitude, (a) ^ A ^ XoA»>o (no noise). 

(b) (xi/«.) < A < (x',/t»o) Vl + x*/4 (primary noise), (c) A ^ (zJ/m.)V1 + »*/4 
(secondary noise). 

occurrence of half-cycles, which is "sbt/bi/ir. The total noise power for 
the two-band process is then the average of this quantity over the 
distribution on A, i.e., 



1 lb 2 



(15) 



<S(A)p(A)dA. 

(pdfof A) 

Since .4 is known to be a Rayleigh variate of mean-square value 
(2bi/b2), the only unknown is the energy function 8(A). 

3.2 Derivation of 8(A) 

Three distinct regions of A can be identified, each giving rise to a 
distinct pattern of signal and noise. The region depicted in Fig. 4a 
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Fig. 5 — Noise energy per half-cycle as a function of A. 



(Region 0) is one in which no slope overload occurs because | x' (t) | 
< x' at all t. Clearly, 6(A) =0 in this case. The region depicted in 
Fig. 4b (Region 1) is one in which slope overload occurs in each half- 
cycle, but over less than the complete interval. We see that the noise 
in this case is primary noise, as denned earlier. Finally, the region 
depicted in Fig. 4c (Region 2) is one in which slope overload occurs 
for the entire duration of each half-cycle. The noise in this case is 
seen to be secondary noise.* 

To find & in Region 1, we analyze the burst spanning [t„, t + At] 
in Fig. 4b. The noise energy in this burst is 



8 = 



t +M 



[A cosw„< — [x — x' (t — t ))~\ 2 dt. 



(16) 



The quantities x„ and t are found by means of trigonometric identities 
to be 



1 



x = VA 2 - (x'o/uo) 2 and t„ = —sin" 1 

03 



\AuJ 



(17) 



Unfortunately, the quantity At cannot be solved for explicitly, but is 



* For the special process under discussion, we see that a given sample function 
has either no noise at all, primary noise onlv, or secondary noise only. For more 
spectrally distributed processes, both kinds of noise occur in the same sample func- 
tion and in a less regular pattern. 
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related to A by the transcendental equation 

Fio At - sin (u„At)~] 2 fAu \ 2 A 

L l-co,(^,) J " {-*) - ' - f - (18) 

We can relate 8 to A (or £) by combining (16), (17), and (18) and 
eliminating the common parameter, w At, between (16) and (18). 
The result is a unique correspondence between ulS/(x'„) 2 and £ having 
the variation shown in Fig. 5 for Region 1. (Note that Region 1 cor- 
responds to the ^-interval [0, tt 2 /4].) A highly accurate functional 
description of this variation has been found to be 



u 8 81 



r /2 [0.43 exp (-0.5014$) + 0.57 exp (-2.10$)]; 



(x' o y 140x 

^ £ < tt 2 /4. (19) 

Comparison with exact results shows this function to be accurate to 
within 0.8 percent. 

To find 8 in Region 2, we analyze the burst spanning [t , t + At] 
in Fig. 4c. In this case, x and t are given by (17), as before, but 
(ji At is precisely tt for all A. Applying these relationships to (16) and 
performing the indicated integration, we obtain 

u„8 1 f / tt 2 \-\ 

This variation is also shown in Fig. 5. 

3 . 3 Expression for N,b (S) 

We now have expressions for 8 in the three regions of A and can 
average this complete result over the distribution on A (or £) to find 
Ntb, (15). By using (18), along with u\ = b 2 /b u S = x' /^b t , and the 
fact that A is a Rayleigh variate of mean-square value (26?/6 2 ), we 
obtain 

, x S * I s *\ I S 2 \ 
p(£) = — exp ( - — exp £ ; £ ^ - 1 



2 \ 2 

= 0; elsewhere. (21) 

Combining this with (15), (19), and (20), we obtain 

N, b(S) - Fm + (|)^) = gj) F( S) (22) 

Primary Secondary 

Noise Noise 
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where F(S) & Fi(8) + F 2 (S), and 



Fi(S) = 
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560 



(i 



-)5*exp(-- 



S- 



0.43Q 



— h 0.5014 

4 \2 



Jx 2 /S 2 
+ 0.57-2 {- 



(f + 210 )|] ; 



F t (S) = 



Q{u\ = 



/5tt 2 -36' 

1 4- ( )S 2 



,9/2 



24 
105 V^ 



exp 



- (1+I 



& 



(23) 
(24) 



16 



erf (Vm) 



[■ 



7 35 

| w 7/2 + - W 5/2 _| w 3/2 + 

2 4 



105 "I 

w ,/2 

8 J 



exp(-u) • (25) 



The variation F(S) is shown in Fig. 6 along with its component parts, 
Fi(S) and F 2 (S). These curves indicate, for the two-band process, 
which region of S is dominated by primary noise and which region 
by secondary noise. 

The only inexactness in our result for N, b (S) arises from the func- 
tional fit to S in Region 1, (19). Since this fit is accurate to within 
±0.8 percent at all £, and both S($) and p(£) are non-negative func- 
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Fig. 7 — Comparison of Nu(S) with previous results. 



tions, the given expression for N t b(S) is also accurate to within ±0.8 
percent (or ±0.035 dB). 

3.4 Comparisons and Interpretations 

The variation of iV,& with S is plotted in Fig. 7, along with the varia- 
tions of N R , (9) ; N P , (10) ; and N A , [(14) with K = 1]. Of particular 
interest is the fact that N t b(S) converges with the results of Rice and 
Protonotarios as S — -> » . The approximations underlying the analyses 
of Rice and Protonotarios become increasingly valid as S increases, 
so the convergence of N t b(S) with their results is not surprising. The 
apparent dependence of noise power on the moments bi and b 2 alone 
at large S is explained by the fact that slope overload bursts are of 
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short duration in this region ; hence, they are shaped primarily by the 
lower-order curvature of x(t), which is reflected in bi and b 2 . As S 
decreases towards small values, however, the higher-order curvature 
of x(t), reflected in b 3 , b h etc., also influences the noise bursts and, 
consequently, the noise power. We should therefore expect the true 
noise power [AT(S)] associated with a given spectrum to reflect other 
features of that spectrum at low S besides 61 and b 2 . If this is so, then 
N t b(S) cannot be assumed to be a universal formula. 

To show that this expectation is correct, we note that, as S — > 0, 
N lb converges to b\/b 2 , which is precisely the ac power (<x 2 ) in the two- 
band process. This is a quite general result, i.e., N(S) always con- 
verges to <r 2 as S — » 0. However, some other process having the same 
b , 61, and b 2 can have an ac power anywhere between bjb% and b , 
so that N(S) will not converge to b\/b 2 for all spectra having these 
moments in common. 

We conclude, then, that N(S) and N lb (S) converge at high S but 
become increasingly dissimilar as S decreases towards zero, the dis- 
similarity depending on the differences between the actual process 
spectrum and the two-band spectrum having the same b , b u and b 2 . 
In Part 2 (following), we will derive a factor which relates N lb (S) 
and N(S) when X(f) is not a two-band spectrum. 

Part 2 

IV. NOISE POWER FOR S ^ 3.5 

We now derive a noise power formula applicable to all spectra, with 
S large (primary noise only). The derivation combines the power series 
method of Rice with the infinite-time averaging procedure used by 
Zetterberg and Protonotarios. The differences here are that many 
terms of the power series expansion are included, and steps are taken 
to minimize the contributions from spurious bursts. The complicated 
expression that results is reduced to a simple formula [#*(£)] that 
displays explicitly the influence of the higher-order spectral moments. 

4.1 Power Series Representation of the Noise Burst 

Let t = be the time origin of a particular positive-going burst, so 
that x'(0) = x' and x"(0) > 0, e.g., Fig. 2. Beginning at t = 0, y{t) 
follows the straight line x(0) + x' t until this ramp function intersects 
x{t) again. We can therefore define the quantity 
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v(t) = x(t) - O(0) + x' Q; t > (26) 

and write the variation of the noise burst as 

n(t) = v(t) if v(t') > for all < t' ^ t 

— otherwise. (27) 

Henceforth, we shall use the symbols n t and v t to denote n(t) and v(t), 
respectively. 

Following Rice, we can express x{t) by a power series expansion 
about t = 0. Noting that .t'(0) = x' and denoting x"(0) by x", v t can 
be written as 

1 1 

vt = hx'W + — x"'(0)t 3 + x^(0)t m + • • • . (28) 

3! ml 

We now invoke the property of Gaussian random processes that the 
?nth and (m — 2)th derivatives at a given time instant are related by 

x (m) _ x ( m -2) _|_ dm (29) 

fem-2 

where d m is a zero-mean Gaussian variate of mean-square value 

d n = b m - (30) 

Om-2 

With this relationship we can write v, as follows : 



v = ) x "i 1 i._ .. 

I °L2! 6 2 4! b 2 b< 6! 

' T— — - — bi t& hi bi b * v 11 

X ° L&i 3! 6^5! biT t 'T 6 7\ J| 

[ t 3 t l rb 4 -i t 5 rb s i t 6 

I 3! 4! L&a J5! Lb 4 J 6! 
- A ( + S t (31) 

where i4, is the first bracketed quantity and 5, is the second. We see 
that A t is the conditional mean of v t given that x'(0) = x' and 
z"(0) = x'o ; and that h t is a linear sum of mutually independent zero- 
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mean Gaussian variates (d 3 , d i} etc.) whose mean-square value at 
time t is therefore 



7; 



♦(Ha--)"*--} - 

We thus have the result that the conditional pdf of v at time * is 

pw^'i—exp |-jp7r)T (33) 

4.2 Derivation of Mean Noise Power 

Let n? be the conditional mean square value of n at time i, given 
that s'(0) = x„ and x"(0) = x' ' >_0. To find the mean burst energy 
when x'(0) = xj, we must average n* with respect to x" and integrate 
the result over all time. We can then multiply this mean energy by the 
mean rate of primary noise burst occurrences to obtain the mean noise 
power. This approach is valid so long as secondary noise can be 
neglected, which it can in the region S ^ 3.5 under consideration. 

From Protonotarios, 3 we know that the correct density function for 
averaging with respect to the conditions x' (0) = x' and x" > is 



p(x' ') = — exp 
o 2 



(*o) 2 



2b, 



x'J ^ 



(34) 



while from Rice 8 we know that the mean rate of noise burst occur- 
rences is (V62/&1/V) exp(-S 2 /2). The noise power in the region S ^ 3.5 
can thus be accurately given by 



1 6 2 / <s 2 \ r 



Xo 

~b~2 
■exp 



(tf)' l 

26 2 I 



nKx' ,x'J)dx' 'dt. (35) 



The remaining analytical task is to find n t {x' , x' '). 

The method of Protonotarios in this regard amounts to averaging v t 
over positive values using the conditional pdf of v t , (33). To do this, 
however, is to incur the spurious burst problem depicted by Fig. 8a. 
The following line of reasoning leads to an improved procedure. 
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\— * l 



(a) 



v,,v,' 




Fig. 8 — An illustration of v t and n t . (a) Vt and n,. (b) v t and yj 4) . 

Let A t denote the partial power series for A t up to and including 
the fourth-order term, and let 



(4) t (4) , . 
v, = A, + 8 t . 



(36) 

We can show that vj 4) ~ v, when t is small and that v t > v ( t 4) in general. 
What we have found in particular is that, for S greater than about 2, 
the difference between v t and v, (4) does not widen appreciably over a 
typical burst duration (see Fig. 8b). Also, since A\ A) is a strongly de- 
creasing function of t beyond its peak, v ( t 4) tends to be the same. 

We therefore calculate nf according to the following criterion: A 
given value of v t is counted as part of the noise (i.e., h, is assumed to 
be v,) if either [v, > 0, v, (4) > 0} or \v t > 0, v { t 4) < 0, v t < 0} ;* other- 
wise, v, is not counted and n, is assumed to be 0. For the sample func- 
tion shown in Fig. 8b, this means that the energies in intervals A and 
B arc both counted in the calculation of »if; the energy in interval C 
is not counted, since y, (4) < and v, > ; and the energy in interval D 
is counted, even though it is clearly spurious. We can reduce contribu- 

* The symbol t), stands for dv(t)/dt. 
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tions of the latter type, however, by using a suitably finite upper limit 
in the integration over time in (35). That is, we choose an upper limit 
sufficiently high that virtually all legitimate noise contributions are 
counted in the calculation but sufficiently low that some spurious 
contributions are omitted. By studying the time variations of A , and 
8 t under a wide variety of conditions, we have settled on an upper 

time limit of 3.5 ^b~Jb 2 . 

Using the above criterion, we have been able to show that 

2 



n 2 = 
t 



— 1(1+ 2a 2 )[l + erf (a 4 )] + -= (2a - a 4 ) exp ( -a!) 

2 I Vtt J 

1 I 

+ - (erfc (d)) (1 + 2a 2 )[erf (a) - erf (o 4 )] 



2 2 

+ — [a exp (-a 2 ) - (2a - a 4 ) exp (-a 4 )] 

V7T 



(37) 



where : 



and 



a = A,/V27 2 ; (38) 

ai = aT/^M; (39) 



To complete the derivation of N(S), we insert this complicated expres- 
sion into (35) and perform th e inte grations over x" and t, remembering 
to use an upper limit of 3.5 V6i/6 2 in the time integral. 

4.3 Simplified Formula 

The formal solution just described is believed to be more exact 
than that of Protonotarios [eq. (53) of Ref. 3], because it entails a 
substantial reduction in the spurious burst contributions. Furthermore, 
in reducing his formal solution to computation, Protonotarios makes 
a number of functional approximations which obscure the influence 
of the higher-order moments b 3 , b h etc. In reducing the present solu- 
tion to computation, we must also make approximations but of a 
different kind, i.e., those associated with performing numerical double 
integration and truncating the infinite power series in the integrand. 
The impact of the higher-order moments, however, is not approxi- 
mated away in the process ; in fact, the following discussion develops 
an empirical expression [denoted by Ni(S)2 which exhibits these 
influences explicitly. 
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We define a dimensionless quantity, which we call the nth-moment 
excess, as 

a b„ - (6n_l/6 n _ 2 ) 

M n = -; n>3. (41) 

2 1 

From (6) we know that M n ^ for all n ^ 3 and, using Fig. 3, we 
can show that M n = for all n ^ 3 if and only if X(f) is the two- 
band spectrum. We can therefore say that N(S) differs from N lb (S) 
only to the extent that M 3 , M*, etc., are not zero. A simple model 
that reflects this fact is one which describes the logarithmic difference 
between TV and A^u as a linear combination of the moment excesses, 

\n(N/N lb ) = AM 3 + 5M 4 + CM 5 + • • • (42) 

where the coefficients A, B, etc., are, in general, functions of S. To 
the extent that this representation is accurate, it is reasonable to 
assume further that the higher-order terms, involving M 6 , M lt etc., 
are negligible for S ^ 3.5. The reason is that these moment excesses 
have little influence on slope overload bursts at such high slope- 
following capacities. The result of this line of reasoning is an approxi- 
mation of the form 

N(S ^ 3.5) = N lb (S) cxp\A(S)M 3 + B(S)M A + C(S)M b ] 

& Ni(S). (43) 

We have tested this approximation by applying the double integral 
solution for N(S), (35) and (37), to a number of spectra, using a 
ninth-order power series to represent v,. With the computed noise 
powers for these spectra, we have then obtained results for j4(&), 
B(S), and C(S) over the range 3.5 ^ S ^ 10.0 (Fig. 9) which appear 
to be quite accurate. In particular, the use of (43) with these results 
is found to predict N(S ^ 3.5), as computed from (35), to within 
0.2 dB for all spectra of practical interest. 

There exist, however, special conditions on X(f) for which (43) 
yields too high estimates of the true noise power. This can occur when 
X{f) contains an isolated high-frequency component which contributes 
materially to 63, b A , and/or b h but not to b , b u or b 2 . Under such circum- 
stances, x(t) might change too rapidly during the noise bursts to be 
properly analyzed by the power series approach as used here. Calcu- 
lations indicate that these cases fall outside the range of practical 
interest, so we need not consider them further. 
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Fig. 9 — Empirical results for A, B, and C; S ^ 3.5. 



V. NOISE POWER FOR S ~ 

We now derive another noise power formula applicable to all spectra, 
but with S very small (secondary noise only). The representation we 
seek is a partial power series of the form 

N(S « 0) = o- 2 [C + CiS + C 2 S 2 )] £ tf.OS) (44) 

where a 2 is the ac power of the input process. 

5.1 Formulation 

The analysis assumes [x(t)\ to be just the ac part of the input 
process, so that x 2 (t) = <r 2 . Since dc components have no effect on noise 
power for S ^ 0, we are justified in ignoring them. 

We consider x'„ to be so small compared to Vbi (i.e., S so small 
compared to unity) that the delta modulator is always in slope over- 
load. In this case, the decoded output [denoted here by y(t) rather 
than y(t)2 is just a succession of alternating ramps of slope ±x„, the 
slope polarity reversing whenever y(t) intersects x(t). This situation is 
illustrated in Fig. 10b. A valid model of the delta modulator for 
analyzing this case is given by Fig. 10a (cf. Ref. 9), from which we 
see that the noise signal is 

n(t) = x(t) - y(t) = x(t) - x' / n q (u)du (45) 



where 



n,(-) ^sgn{n(-)}. 



(46) 
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Fig. 10 — Equivalent delta modulator and signals for S « 0. (a) AM model for 
S ~ 0. (b) Input and feedback signals. 

The noise power in the region S = can thus be accurately given by 



N(S) = x*(t) - 2x' 



x(t) I n q (u)du 



-2x(t)y(t) 



+ (x'oV 



n q (u)n q {v)dudv. (47) 



y\t) 

If we envision each of these three terms as a power series in x' , we 
see that x 2 (t) has a zero-order term only (i.e., is independent of x' ) ; 
— 2x(t)y(t) may contain first-, second-, and higher-order terms; and 
y 2 (t) may contain second-order terms and higher. We can thus rewrite 
(47) as 

N(S) = K + ZKiX', + K 2a (x' ) 2 +•••] + lK 2b (x' Y +•••]• (48) 



xKt) -2x(t)y(t) y*(t) 

Retaining only terms up to second-order, substituting x' = Svoi, and 
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comparing with (44), we see that 

a 2 a 2 <t 2 

Clearly, K = x 2 (t) = a 2 , so that Co = 1 in all cases. The remainder 
of this section is devoted to finding C\ and C 2 . 

5.2 Derivation of Ci 



We will derive C\ by evaluating the first-order term of — 2x(t)y(t), 
(47), and then invoking (49). For x' very small, n q {t) does not differ 
appreciably from x q {t) = sgn(.r(0 ), which is the quantizer response to 
x(t) applied alone. It is therefore fruitful to represent n q {t) as 

n q (t) = x q (t) + 9(t) (50) 

where 9(t) is a correction signal related to the finiteness of x' and is 

defined by 

r+2 if y < x < 

d = \-2 if y > x > (51) 

I otherwise . 

The relationships between n q (t), x q (t), and 0(t) are illustrated in Fig. 
11. Although the peak magnitude of 6(t) is fixed, its "duty cycle" can 
be seen to depend directly on x' . 

Inserting (50) into the cross term of (47) and interchanging the order 
of integration and averaging yields 

-2x{t)y(t) = - 2x' / x(t)x q (u)du - 2x' / x(t)6(u)du. (52) 

J —oo J — « 

Since the duty cycle of 6{t) vanishes as x' -> 0, we can see that the 
first term must be identical to Kix' in (48), while the second term 
contains K 2a (x' ) 2 plus higher-order terms in x' . 

To evaluate the first term, we use a relationship applicable to 
[x(t)\ because it is a Gaussian process, namely, 



x(t)x q (u) = -v/ = \ (53) 

where R x (t) is the autocorrelation function of \x(t)\. We now apply 
this to the first term of (52) and use the fact that the area under 
R z {r) from t = to t = °° is jX(0). Equating the result to K\x' 
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Fig. 11— The relationships between x(l), y(t), x„(t), n q {t), and 0(0. 



we obtain 



which, from (49), yields 



K,= 



1 X(0) 



i rx(0)V6fi 



(54) 



(55) 



5.3 Derivation of C 2 

The general approach used in the preceding analysis was also applied 
to finding C 2 . Unfortunately, it leads to an approximation which is 
both highly complicated and not sufficiently precise for many spectra. 
For this reason, we have resorted to finding C 2 for a not-quite-Gaussian 
process having the same spectrum as the specified Gaussian one, and 
verifying that it is virtually exact under what appear to be "worst- 



412 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1973 

case" conditions. The details, which are documented but not pub- 
lished, 10 are quite involved and so we shall merely outline the approach. 
Consider a wide sense stationary process \x(t) } having sample func- 
tions of the form 

x(t) = A cos (art + 0) (56) 

where A is a Rayleigh-distributed amplitude of mean-square value 
2<r 2 ; is a random phase uniformly distributed over [ — ir, +7r]; and 
w is a random frequency whose pdf is 

p(C) = X(f = |«|/2ir); all*. (57) 

It can be shown that the power spectrum for this process is X(f), and 
that the ensemble pdf of x at any t is Gaussian with zero mean and 
variance a 2 . Thus, \x(t)} has certain properties in common with the 
Gaussian process \x(t)} having the spectrum X(f). 

We can derive the mean slope overload noise power for the process 
[x(t)\ by finding the noise power for a single sample function, (56), 
and averaging over the distributions on A, </>, and w. Using the methods 
of Section III, we can show that the averaging over A and <t> yields 

(y/2x' \ 
— - (58) 

VA 2 aV 

where F(-) is defined by (23)_through (25). Averaging (58) over the 
pdf of a, (57), and replacing A 2 by 2<r 2 and x' by S^bl, we obtain 

ft(8) = rX(f)F(-^)df. (59) 

The derivatives of this function with respect to S can be found quite 
simply, although some caution is required when evaluating them at 
S = 0. The resulting value for C 2 ( = £#"(0)) is found to be 

itrrotw^/Mj (60) 

s-,0 \Sa 3 J 2irf dS \2irfaJ j 

For most realistic processes, the continuous part of the power spec- 
trum has zero slope at / = 0, i.e., X'(0) = 0. In that circumstance, 
(60) reduces to the much simpler form 

\12 /a* J* (2tt/) 2 
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Fig. 12 — Two families of spectra, (a) Truncated Butterworth spectrum, (b) Band- 
pass uniform spectrum. 



In general, the above C 2 for the non-Gaussian process {55(0) can- 
not be compared with the C 2 for the Gaussian process \x(t) \ since we 
do not have a result for the latter case that is good for all X(f). We 
do, however, have exact results under some special conditions on X(f) 
for which C 2 is expected to be maximally dissimilar for the two pro- 
cesses. The very close agreement observed under these conditions 
persuades us that (60) is a reliable representation of C 2 for Gaussian 
processes having any A'(/). 

5.4 Final Result for N,(S) 

Our final approximation to N(S ~ 0) is (44) with Co = 1 and d 
and C 2 given by (55) and (60). Note that when X(0) = 0, Ci is zero 
and C 2 is negative, so that N,(S) is convex at the origin; and that when 
X(Q) ?* 0, C\ is negative and C 2 can be either positive or negative, 
depending on X(f). For a spectrum like the one in Fig. 12b, this 
implies a functional discontinuity at /3 2 = 0, which can be explained 
as follows: So long as 2 ^ 0, the curvature of N S (S) is convex at the 
origin, but it becomes concave at some S > 0. As 2 decreases, the 
curvature at the origin becomes sharper and the point of inflection 
occurs at successively lower values of S. Finally, in the limiting case 
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02 = 0, the point of inflection occurs at S = and so N, (S) becomes 
concave at the origin. Thus, an abrupt increase in X at / close to 0, 
reflected in a large negative value of C 2 , signifies a very small range 
of validity and calls for caution in using the new result for N(S « 0). 
For concreteness, suppose that the continuous part of X(f) in- 
creases to a large peak density at some low frequency f , where 
jo « V&i/27T(r. To obtain a more useful approximation to N(S ~ 0) 
under such circumstances, Ci and C 2 should be recomputed as if the 
point of high density were at / = instead of / = /„, i.e., as if the 
spectrum were X(f + /„). The resulting N,(S), (44), then corresponds 
to a spectrum quite similar to X(f) but without the abrupt change 
at low /. Consequently, it exhibits the correct "large-scale" curvature 
at low S while omitting the sharp "small-scale" curvature at S near 
associated with the true C\ and C 2 . An empirically derived rule-of- 
thumb that leads to good results in all cases is the following: If 
C 2 < - 10, shift the continuous part of X(f) to the left until the 
nearest high-density peak is at / = 0, and recompute C\ and C 2 for this 
modified spectrum. Then compute the quantity C = (0.5 C 2 — 10 Ci) 
using these values. If C exceeds the original value of C 2 , use the newly 
computed values of d and C 2 ; otherwise use the original ones. 

VI. GENERAL NOISE POWER FORMULA 

We now seek an expression for slope overload noise power [N <>(£)] 
which is accurate over all S for all spectra and is simple to use. Our 
approach is to derive separate functional descriptions for the regions 
< S < 4.0 and S ^ 4.0. The descriptions are such that the resulting 
N (S) and its first derivative are continuous at the boundary S = 4.0. 

6.1 N (S)forS ^ 4.0 

The result Ni(S), (43), gives a very accurate approximation to N(S) 
in the region S ^ 3.5. However, our result for N tb (S), (22) through 
(25), is quite complicated, and the results for A (S), B(S), and C(S) 
are available in graphical form only (Fig. 9). We now use these results 
to derive a simple but still accurate representation for this region. 

Taking a suggestion from Rice's analysis, 2 we speculate that N(S) 
for S ^ 4.0 can be accurately approximated by 

243 6 2 ie xp(->SV2) 
N(S ^ 4.0) = — 7= — 

•[exp (-a ie xp (-a 2 £))] = N (S ^ 4.0). (62) 
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The quantity preceding the brackets is Rice's result N R (S), (9), while 
the bracketed term (with a 2 > 0) is a monotonically decreasing func- 
tion that converges to unity as S — > ». We choose ct\ and a 2 so that 

tf .(4.0) = iV,(4.0) and iV' o (4.0) = ATJ(4.0) . 

By using the data in Fig. 9 to obtain A (S), B(S), COS), and their first 
derivatives at S = 4.0, and the data in Fig. 6 to obtain N t b(S) and 
its first derivative at S = 4.0, we obtain 

0.2141 - 0.024il/ 3 - 0.196M 4 - 0.067 ilf 5 

«2 = — ■ (63) 

0.6394 - 0.057M 3 - 0.426-/tf 4 - 0.093M 5 

and 

ai = (0.6394 - 0.05771/ 3 - 0.426M 4 - 0.093ilf B ) exp (4a 2 ) (64) 

where M%, M 4 , M b are the moment excesses denned by (41). Note that 
when X(f) is the two-band spectrum, (M„ = 0, n ^ 3), on and a 2 
reduce to 2.45 and 0.336, respectively. 

The above approximation to N(S ^ 4.0) can be tested by comparing 
it with the more precise results computed for various spectra using 
(35) and (37). The indications from such comparisons are that the 
approximation is accurate to within ±0.3 dB over S ^ 4.0 for all 
spectra. 

6.2 N (S) forO < S < 4.0 

The analyses of Sections IV and V inform us about N(S) at the 
extremities of the region < S < 4.0, but not about the variation in 
between. To estimate this variation reliably, it is convenient to express 
N(S) in the form 

N(S) = G(S)N tb (S) (65) 

and then seek an approximation to G(S). The latter is a spectrum- 
dependent function that differs from unity because-and to the extent 
that-A'X/) differs from a two-band spectrum having the moments 6 , 
b h and b 2 (Fig. 3). From physical reasoning given earlier, we expect 
that G'(*S) — ► 1 as S — > <x> and that, as S decreases from high values 
towards 0, G(S) increases because-and to the extent that-6 3 , b 4 , etc., 
exceed the minimum values, (6), corresponding to the two-band spec- 
trum. The maximum value of G(S) should then be its value as S — * 0, 
which is u 2 {b 2 /b\). 

This reasoning is supported by careful scrutiny of the results for 
Ni(S) and N,(S), Sections IV and V. From Section IV we can show 
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that G(S) is a nonincrcasing function* for 5 = 3.5 and that it goes to 
unity as 5 — > « . From the results of Section V we can show that (7(5) 
is a nonincreasing function* at 5 ~ and that it goes to a (62/&1) as 
g _» 0. A logical consequence of these observations is that, for < S 
< 3.5, (7(5) is either a nonincreasing function, or has at least two 
extrema. The latter possibility has no physical basis in fact, and so we 
conclude that G(S) is a nonincreasing function over all 5, approaching a 
maximum value of a 2 {b 2 /b\) as 5 — ► and a minimum value of 1 as 
5-» ». 

Computations show that the decrease in G(S) over < 5 < 4.0 is 
less than 10 dB for all spectra of possible practical interest. Thus if we 
can find a functional approximation to the quantity 



a. ..._. , [N(S) 
= In - 



g(S) = In (7(5) 



LN lb (S). 



< 5 < 4.0 (66) 



that is accurate to within ±10 percent, the resulting approximation 
to 0(8) will be accurate to within ±1 dB for all spectra of interest. 
This accuracy should be possible to achieve since g(S) near 5 = 
and 5 = 4.0 are known from the results of Sections III, IV, and V, 
and the above argument persuades us that g(S) is nonincreasing over 
the region in between. 

The function that we use to approximate g (S) is 

g (S) = In (ab 2 /bl) + <nS + a 2 [exp (a 3 S + a 4 £ 2 ) - 1] (67) 

where Oi, • • • a 4 are chosen to give ^(0), ^'(0), g (4), and p£(4) the 
values predicted by the results of Sections III, IV, and V. These values 
are found by using (66), with N. (S) and Ni(S) replacing N(S) at 
5 = and 5 = 4, respectively. Thus, the following equations must be 

satisfied : 

eti + 0,20,3 = Ci\ (68) 



a 2 (a 2 3 + 2a 4 ) = 2C 2 - c\ + f 4 - M; 

In (<r\/bi) + 4ai + a2[exp (4a 3 + 16a 4 ) - 1] 



and 



(69) 
= 0.057M 3 + 0.426M" 4 + 0.093M 5 ; (70) 



ai + a 2 (a 3 + 8a 4 ) exp(4a 3 + 16a 4 ) 

= - 0.024M 3 - 0.196M 4 - 0.067Mb (71) 



* In fact, G(S) is a decreasing function for all but the two-band spectrum, for 
which case it is precisely 1. 
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Table I — Formula Coefficients for the Butterworth Spectrum 



01 


ai 


«2 


a 3 


a 4 


0.02 


-0.177 


2.26 


-3.93 


-8.76 


0.068 


-0.151 


1.29 


-3.70 


-7.04 


0.10 


-0.137 


1.06 


-3.75 


-6.99 


0.25 


-0.096 


0.66 


-3.90 


-7.03 


0.50 


-0.066 


0.50 


-3.91 


-6.67 


1.0 


-0.047 


0.42 


-3.86 


-6.22 


2.0 


-0.040 


0.38 


-3.84 


-6.00 


00 


-0.036 


0.37 


-3.83 


-5.90 



where C,, C 2 and M„(n ^ 3) are given by (55), (60), and (41), 
respectively. 

A set of solutions that is quite accurate in most cases is the following : 



ai = - 0.024M 3 - 0.196M 4 - 0.067M 5 

a 2 = In (a 2 b 2 /bl) - 0.155M 3 - 1.210M 4 - 0.361M 5 
r(Ci — ai)/a 2 ; a 2 > 

a 3 = -j 0; a 2 = 

L|Ci — ai|/a 2 ; a 2 < 



a>\ — i 



2C 2 - C x + ( 4 

6 



a.. 



2 

- a 3 



a-t > 0. 



(72) 
(73) 

(74) 



(75) 



0; a 2 ^ 

The results for a 3 and a 4 when a 2 ^ are approximations only, but a 2 
is nonpositive only when X(f) is relatively narrowband, i.e., when 
<r 2 b 2 /bl m 1, in which case G(S) is close to 1 over all S. Hence, these 
approximations lead to a quite accurate representation of g(S). In the 

Table II — Formula Coefficients for the Uniform Spectrum 



Pi 


Oi 


at 


a 3 


a 4 





-0.036 


0.37 


-3.83 


-5.90 


0.05 


-0.036 


0.32 


-4.44 


-8.21 


0.10 


-0.036 


0.27 


0.14 


-11.96 


0.20 


-0.035 


0.17 


0.21 


-7.58 


0.30 


-0.033 


0.08 


0.41 


-7.95 


0.40 


-0.029 


0.02 


1.44 


-18.30 


0.50 


-0.024 


-0.01 


-1.76 





0.80 


-0.004 


-0.01 


-0.43 





1.0 
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Fig. 13 — Noise power results for several spectra. 

more usual case where a 2 > 0, the above solutions are highly accurate 
so long as 4a 3 + 16a 4 ^ — (4 + lna 2 ), a condition satisfied for most 
spectra. If this inequality is violated, solutions for a h ■ • • a 4 can be 
found by graphical or computerized methods, but such situations 
appear to be rare. 

With g(S) approximated as above, N„(S) over < S < 4.0 can be 
expressed as the product exp (g (S))N lb (S). To obtain a usable ex- 
pression, however, N lb (S) should be given in a more simple form than 
eqs. (22) through (25). An approximation to N, b (S) which is accurate 
to within ±0.3 dB is 



6? 
Ntb(S) = — [1 + 2.753S + 2.952S 2 ] 

•exp (-2.753S - 0.34 IS 2 ); 



< S < 4.0. (76) 
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Fig. 14 — Noise power results in greater detail. 

Combining (76) with (67) leads to the following : 

No(S) = <r 2 [l + 2.7535 + 2.952S 2 ]exp(-0.341£ 2 ) 

•exp { ( ai - 2.753)5 + a 2 [exp (a 3 S + a<S 2 ) - 1]} ; 

< S < 4.0. (77) 
6.3 Final Expression 

A highly accurate two-region approximation to N{S) is given by the 
combination of (62) and (77). We have determined, however, that 
(77) alone predicts N(S) with an accuracy of 1 dB or better for S 
between 4.0 and 6.5 for all spectra. Moreover, we have determined that 
N(S) is at least 119 dB below a 2 for S ^ 6.5 and for all spectra. Ob- 
viously, then, (77) can be used over all S for which N(S) is not negli- 
gibly small. We therefore present (77) as our final expression for slope 
overload noise power, with a h • • • a 4 defined by (6S) through (71) 
and accurately approximated by (72) through (75) in virtually all 
cases. 



VII. NUMERICAL RESULTS 

The new result has been applied to the two families of spectra shown 
in Fig. 12. The Butterworth spectrum is characterized by an upper 
truncation frequency (W) and the ratio (0,) of the 3-dB corner fre- 
quency to W. Note that, when ft — ► oo , the spectrum approaches that 
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Table III— N p /N in dB for the Butter worth Spectrum 



X 


0.036 

(0, = 0.02) 


0.111 
03 1 = 0.068) 


0.289 
(0i = 0.25) 


0.556 

(0i = °°) 


1.0 
2.0 
3.0 
4.0 
5.0 
6.0 


+2.85 
+0.41 
+0.43 
+0.96 
+ 1.14 
+0.90 


+ 1.76 
-0.20 
-0.01 
+0.74 
+ 1.12 
+0.85 


+ 1.44 
-0.22 
-0.04 
+0.71 
+ 1.11 
+0.80 


+ 1.70 
-0.10 
+0.22 
+0.81 
+ 1.13 
+0.70 



of bandlimited white noise. The bandpass uniform spectrum is charac- 
terized by an upper truncation frequency (W) and the ratio (0 2 ) of 
the lower truncation frequency to W. Note that p 2 = corresponds to 
bandlimited white noise in this case and that, as 0* — ► 1, the spectrum 
becomes infinitely narrowband. 

For each of these two families, the formula coefficients, (oi, a 2 , a 3 , a 4 ), 
have been computed as functions of the respective /3-parameter. The 
results for the Butterworth spectrum are given by Table I and those 
for the uniform spectrum by Table II. The values shown are rounded 
to the decimal accuracy required. 

Curves of N /<x 2 versus S for several Butterworth spectra of practical 
interest and for the narrowband spectrum are shown in Fig. 13. (The 
curve for the narrowband case represents an upper bound on N(S)/a 2 
for all spectra.) Also, an informative closeup that treats more spectra 
but over a reduced range of S is given by Fig. 14. 

Since the new formula is highly accurate for all spectra, it can be 
used to estimate the accuracy of the Protonotarios formula, (10). 
Table III compares N P (S) with N (S) for various Butterworth spectra 
spanning the range of y = b\/b 2 b Q of practical interest, while Table 
IV gives the comparison for various two-band spectra. The larger 
errors in Table IV must be regarded as untypical, since the two-band 
spectrum itself is rather artificial. For more typical spectra, we can 

Table IV—N p /N in dB for the Two-Band Spectrum 



X 


0.111 


0.289 


0.556 


1 


1.0 
2.0 
3.0 
4.0 
5.0 
6.0 


+5.04 
+2.45 
+ 1.97 
+2.05 
+ 1.79 
+0.84 


+3.55 
+ 1.46 
+ 1.23 
+ 1.54 
+ 1.55 
+0.88 


+2.49 
+0.75 
+0.70 
+ 1.12 
+ 1.30 
+0.74 


+ 1.52 
+0.06 
+0.19 
+0.71 
+ 1.02 
+0.50 
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expect the errors in the Protonotarios formula for given <S and 7 to be 
more like the entries in Table III. 

VIII. CONCLUSION 

Aside from providing a new formula for slope overload noise power, 
this work has served to resolve the disparities between previously 
published formulas and can be used to identify the accuracies and 
limitations of those most widely used. Furthermore, the new result 
and the underlying methods of analysis can be extended to treat other 
important topics that have been mostly ignored before. Among these 
are the average duty cycle of slope overload, the slope overload noise 
spectrum, the effects of leaky feedback integrators on slope overload 
noise power, and the slope overload noise power for certain non- 
Gaussian processes. 
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